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Abstract 

We further study the validity of the Monte Carlo Hamiltonian method. The 
advantage of the method, in comparison with the standard Monte Carlo Lagrangian 
approach, is its capability to study the excited states. We consider two quantum 
mechanical models: a symmetric one V{x) = |a:|/2; and an asymmetric one V{x) = 
oo, for X < 0 and V{x) = x, for x > 0. The results for the spectrum, wave functions 
and thermodynamical observables are in agreement with the analytical or Runge- 
Kutta calculations. 
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1 Introduction 


Path integral quantization in the Lagrangian formulation and canonical quantization in 
the Hamiltonian formulation are two ways to quantize classical systems. Either one has 
advantages and disadvantages. The Lagrangian formulation is suitable for numerical simu¬ 
lations on a computer via Monte Carlo (MC). The enormous success of lattice gauge theory 
over the last two and half decades is certainly due to the fact that the MC method with 
importance sampling[1] is an excellent technique to compute high dimensional (and even 
“inhnite” dimensional) integrals. 

Unfortunately, using the Lagrangian formulation it is difficult to estimate wave func¬ 
tions and the spectrum of excited states. Wave functions in conjunction with the energy 
spectrum contain more physical information than the energy spectrum alone. Although 
lattice QCD simulations in the Lagrangian formulation give good estimates of the hadron 
masses, one is yet far from a comprehensive understanding of hadrons. Let us take as 
an example a new type of hadrons made of gluons, the so-called glueballs. Lattice QCD 
calculations[2] predict the mass of the lightest glueball with quantum number = 0’*'^, 
to be 1650 ± lOOMeU. Experimentally, there are at least two candidates: /o(1500) and 
/j(1710). The investigation of the glueball production and decays can certainly provide 
additional important information for experimental determination of a glueball. Therefore, 
it is important to be able to compute the glueball wave function. 

In the Hamiltonian formulation, one can obtain the ground state energy, but also wave 
functions and the spectrum of excited states. Often, and in particular in the case of many- 
body systems, it is difficult to solve the stationary Schrodinger equation. Recently, we have 
suggested how to construct an effective Hamiltonian via MC[3] to describe the low energy 
spectrum and wave functions. The method has been tested in quantum mechanics[3, 4] in 
1,2 and 3 spatial dimensions (the free system, the harmonic oscillator, the quartic potential, 
the local potential V{x) = —Vosech^(x)) and the scalar model in quantum field theory[5] 
(Klein-Gordon model and 0"^ scalar model). 

The purpose of this paper is to further test the viability of the MC Hamiltonian method 
in case of two potentials which generate a bound state spectrum: 


and 


V{x) = \x\/2, 



a; < 0 
a; > 0. 


( 1 ) 

( 2 ) 


These two models have their own interest in quantum theory. They have the form of a static 
quark-antiquark conhning potential occurring in QCD. The second potential describes the 
particle moving above the surface of the ground. In our previous investigations, only 
symmetric potentials were considered. Here we would like to see if the method works for 
linear potentials and those with infinity and asymmetry. 
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The remainder of the paper is organized as follows. In Sect. 2, the basic ideas of the 
effective Hamiltonian method is briefly discussed. In Sect. 3, the numerical results are 
given. The conclusions are summarized in Sect. 4. 


2 Effective Hamiltonian 


The construction of an effective Hamiltonian starting from a regular basis in position space 
has been proposed in Ref. [3]. Following Feynman’s path integral formulation [6], we 
consider the transition amplitude in imaginary time between time t = 0 and t = T. Us¬ 
ing imaginary time makes the path integral mathematically well defined, and renders it 
amenable to numerical simulations by MC. Because the effective Hamiltonian is time inde¬ 
pendent, its construction in imaginary time should give the same result as in real time. We 
consider the transition amplitude for all combinations of positions Xi,Xj G {xi,... ,xn}, 


Mij{T) =< Xi\e >= J[dx]exp[—SE[x]/h] 


Xi,T 


(3) 


where Se denotes the Euclidean action. From the transition amplitudes Mij{T), one can 
construct the matrix 

M(T) = [M,,,{T)]nxn. (4) 

M{T) is a positive, Hermitian and symmetric matrix. M can be factorized into a unitary 
matrix U and a real diagonal matrix D{T), such that 


M(T) = U^D{T)U. 


(5) 


Then from Eq.{3), Eq.{5) one can identify 


UL=<x,\Ef>>, D 


,t(T) = 


T/n 


( 6 ) 


The k — th eigenvector \Elf'^ > can be identified with the k — th column of matrix W. The 
energy eigenvalues E’^^^ are obtained from the logarithm of the diagonal matrix elements 
of D{T). This yields an effective Hamiltonian, 


Heff = 


N 

E 

k=l 


El^^ > EV^ < E. 


eff 


W/| 


(7) 


We compute the matrix element Mij{T) directly from the action via MC with importance 
sampling. For details see Ref. [3]. 

To get the correct scale for the spectrum, the position state \xk) at initial time or 
hnal time tf (Bargman states or box states) should be properly normalized. We denote 
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a normalized basis of Hilbert states as |efc), k = In position space, it can be 

expressed as 


, . . l/VAXfc, xe [Xk,Xk + l] 

= { 0. ® 

where Axk = Xk+i — Xk- There is some arbitrariness in choosing a basis for the initial and 
hnal states. The simplest choice is a basis with = const., which is called the “regular 
basis”. 


3 Numerical results 


For the potential in Eq. (1) an exact solution of the ground state energy and wave function 
is analytically known [7]. They are given by 


rpExact 

-^0 

'4)q{x) 



3/2 


0.5094(hVm)^/^ 


(9) 


1 

where ^ = (|a;| -2Eq) (m/Ti^^ , = ^U=o, |^o| = zo is the solution to the 

equation 


■h/Z^z) — J- 2 / 3 {z) — 0, (10) 

and Ji /3 and J_i /3 are Airy functions. We compare two methods: the MC Hamiltonian 
method and the Runge-Kutta algorithm with the node theorem. The latter one means that 
we solve the stationary Schrodinger equation expressed as differential equation in position 
space and search for the discrete eigenvalue i7„ corresponding to the number of nodes of 
the wave function. Here n denotes the quantum number, being related Eq. (7) to k by 
n = k — 1. Because the transition matrix is symmetric, we implement in the code 

Tfij = {Mij + Mj^i) /2, 

= Mij, (11) 

when calculating the matrix elements using the MC Hamiltonian method. There is an 
additional symmetry for the potential in Eq. (1): V{—x) = V{x) so that the matrix 
elements should satisfy MAr+i_i,Ar+i_j = Mjj. Therefore, in the code 

Mij = {Mij + MN+i-i,N+i-j)/‘2, 

Tf7v+i-j,v+i-/ = Afjj. (12) 


3 



n 

j^Exact 

En 




0 

0.509 

0.509 

0.504 ± 0.002 

1.0% 

1 


1.169 

1.171 ± 0.002 

0.2% 

2 


1.624 

1.622 ± 0.002 

0.1% 

3 


2.044 

2.047 ± 0.003 

0.2% 

4 


2.410 

2.407 ± 0.004 

0.1% 

5 


2.760 

2.764 ± 0.007 

0.1% 

6 


3.082 

3.075 ± 0.021 

0.2% 

7 


3.393 

3.399 ± 0.034 

0.2% 

8 


3.686 

3.669 ± 0.078 

0.5% 

9 


3.972 

3.982 ± 0.019 

0.3% 

10 


4.244 

4.202 ± 0.105 

1.0% 

11 


4.511 

4.532 ± 0.041 

0.5% 

12 


4.768 

4.683 ± 0.069 

1.8% 

13 


5.020 

5.074 ± 0.142 

1.0% 

14 


5.264 

5.155 ± 0.153 

2.1% 

15 


5.504 

5.673 ± 0.137 

3.1% 

16 


5.738 

5.708 ± 0.117 

0.5% 


Table 1: Spectrum corresponding to the potential in Eq. (1). Here and E^'^' 

denote respectively the data from the analytic solution, Runge-Kutta algorithm and MC 
Hamiltonian. The input parameters are m = 1, h = 1, Ax = 1 and iV = 21. The statistical 
an relative errors of E^^ are also shown. The error of the Runge-Kutta algorithm is less 
than the last digit and not shown. 


The use of these symmetries lead to reduction of the numerical errors. We have computed 
the energy spectrum, wave functions as well as thermodynamical observables such as av¬ 
erage energy, specihc heat and the partition function. The results can be summarized as 
follows. Fig. 1 shows the ground state wave function. One observes good agreement be¬ 
tween the MC Hamiltonian result and the exact solution. Fig. 2 shows a similar comparison 
for the wave function of the hrst excited state. Tab. 1 shows the data for the spectrum. 

The model described by the potential in Eq. (2) is exactly solvable. The spectrum and 
wave function are given by 


E„ 




A„(hV2m)^/^, 



3 


0 , 


\e'^) + J-i 



X < 0 

X >0, 


(13) 
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n 

j^Exact 

En 

E^^ 


0 

1.856 

1.774 ± 0.086 

4.4 % 

1 

3.244 

3.182 ± 0.112 

1.9 % 

2 

4.382 

4.346 ± 0.125 

0.8 % 

3 

5.383 

5.300 ± 0.082 

1.6 % 


Table 2: Spectrum corresponding to the potential in Eq. (2). The input parameters are 
m = 1, h = 1, Ax = 1 and N = 20. 


1 /3 

where = {x — En) (2m , and is the solution to the equation 

(14) 

Here are the first few: Aq = 2.338, Ai = 4.088, A 2 = 5.521 and A 3 = 6.787. Fig. 3, Fig. 
4 and Tab.2 show a comparison of the spectrum and wave functions corresponding to the 
potential in Eq. (2). Although the agreement is good, the relative errors are bigger, in 
comparison with the first case. This might be due to the lack of parity symmetry in the 
potential Eq. (2). 

The partition function is dehned by Z{P) = Tr{e~^^). The average energy is dehned 
by U{/3) = Tr{He~^^)/Z, and the specihc heat by C{/3) = dU/dr. We have used the 
following notation: [3 = r is the temperature, ks is the Boltzmann constant, and 

we identify (3 with the imaginary time T hj j3 = T/h. Since we have approximated H by 
Tfeffi we can express those thermodynamical observables via the eigenvalues of the effective 
Hamiltonian 


Z^^{/3) 

U^^{(3) 

C^^{(3) 



(15) 


The numerical results are shown in Fig.5 to Fig. 8 . One observes good agreement between 
the results from the MC Hamiltonian and the reference solution. 


4 Conclusion 

In the preceding sections, we have applied the MC Hamiltonian method[3] to the linear 
potentials: one is symmetric, and another is asymmetric and has inhnity. The resnlts for 
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the spectra, wave functions and thermodynamical observables are in good agreement with 
other methods. We also observe that the use of the parity symmetry whenever exists, can 
lead to systematic reduction of numerical errors. 

In most cases, where no exact solution is available, numerical solution is required. 
Runge-Kutta method works only in ID and it is very difficult to obtain information on the 
excited states by standard Lagrangian MC approach. In our opinion, the MC Hamiltonian 
has an advantage in these cases. 

The MC Hamiltonian method approximates a quantum system with infinite states by 
a hnite system with dimension N. Increasing N, we can obtain more excited states, and 
reduce the systematic errors. Larger N value is also necessary when approaching the smaller 
(3 (or higher temperature regime). However, the computation of the matrix elements is very 
time consuming. In this paper, we used the regular basis, which is good enough for few 
body quantum mechanics and might not be very efficient for quantum held theory. In 
Ref. [5], we suggested the concept of MC Hamiltonian using the stochastic basis and 
applied it successfully to the Klein-Gordon model and 0^ scalar model. We believe that 
the application to more complicated systems will be very interesting. 
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Figure 1: Ground state wave function '^ 0 ( 2 ^) corresponding to the potential in Eq. (1), 
where the continuous line represents the result from the Runge-Kutta method and the 


symbols stand for the data from the MC Hamiltonian method. 
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V{x) = \x\/2 

m = 1,T = l,h = 1, Ax = 1, iV = 21 



Figure 2; Same as Fig. 1, but for first excited state 'ipi{x). 


V{x) = X for x > 0, = cx) for x < 0 

m = 1, T = 1, ^ = 1, Ax = 1, N = 20 



X 


Figure 3: Ground state wave function 'ipo{x) corresponding to the potential in Eq. (2), 
where the continuous line represents the exact solution and the symbols stand for the data 
from the MC Hamiltonian method. 
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V{x) = X for x > 0, K(a;) = cx) for x < 0 
m = 1, T = 1, ^ = 1, Ax = 1, A = 20 



Figure 4: Same as Fig. 3, but for first excited state '^i(x). 

V{x) = |x|/2 

m = l,T= 1,^ = 1, Ax = 1, A = 21 



Figure 5: Average energy U{j3) corresponding to the potential in Eq. (1), where the 
continuous line represents the result obtained by substituting the eigenvalues of the Runge- 
Katta method into Eq. (15) and the symbols stand for the data from the MC Hamiltonian 
method. 
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V{x) = \x\/2 

m = l,T= 1,^ = 1, Ax = 1, A = 21 



Figure 6: Same as Fig. 4, but for specific heat C(/9). 


V{x) = X for X > 0, V^(x) = CX) for X < 0 
m = 1, T = 1, h = 1, Ax = 1, A = 20 



Figure 7: Average energy U{j3) corresponding to the potential in Eq. (2), where the 
continuous line represents the result obtained by substituting the eigenvalues of the anlytic 
solution into Eq. (15) and the symbols stand for the data from the MC Hamiltonian 
method. 
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